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Entanglement entropy of random quantum critical points in one dimension 
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For quantum critical spin chains without disorder, it is known that the entanglement of a segment 
of A'' ^ 1 spins with the remainder is logarithmic in A'^ with a prefactor fixed by the central charge 
of the associated conformal field theory. We show that for a class of strongly random quantum 
spin chains, the same logarithmic scaling holds for mean entanglement at criticality and defines a 
critical entropy equivalent to central charge in the pure case. This effective central charge is obtained 
for Heisenberg, XX, and quantum Ising chains using an analytic real-space renormalization group 
approach believed to be asymptotically exact. For these random chains, the effective universal 
central charge is characteristic of a universality class and is consistent with a c-theorem. 
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Second-order phase transitions at zero temperature 
show universal scaling behavior determined by the col- 
lective physics of quantum fluctuations. Recently the 
scaling of the entanglement near such quantum critical 
points has been of special interest: at a quantum criti- 
cal point, the length scale over which different regions of 
the system are entangled becomes divergent. The entan- 
glement near criticality was shown to obey a universal 
scaling law in some one-dimensional (ID) systems. Most 
quantum phase transitions in pure ID systems are in- 
variant under local conformal transformations, and the 
entanglement at a critical point is related to the central 
charge of the associated conformal field theory [1, 2]. 

Our primary result is that there exists universal entan- 
glement scaling even for a class of disordered quantum 
critical points in one dimension that are not conformally 
invariant. The specific theories considered here describe 
random quantum spin chains: the Heisenberg, XX, and 
quantum Ising chains with random nearest-neighbor cou- 
pling have been previously found, [3, 4] using a real-space 
renormalization-group (RG) approach, to be described 
by strongly disordered critical points, as reviewed below. 

The entanglement of a pure quantum-mechanical state 
1^) with respect to a partition into two subsystems A and 
B is the von Neumann entropy of the reduced density 
matrix for either subsystem: 



-Tip A log2 PA = -Trps log2 pb 



(1) 



where the reduced density matrix pA for subsystem A is 
obtained by tracing over a basis 0^ of subsystem B 

{^a\pa\<P\) = E«<^aI X (<^^IW(^(|0i) X |<^^)). (2) 



Note that this pure-state entanglement of a spin chain [5] 
differs from the two-spin mixed-state entanglement [6, 7], 
which also has special behavior near a phase transition, 
but is only nonzero at short distances and is tied to the 
spin-spin correlator rather than the central charge [8] . 

For conformally invariant critical theories in one di- 
mension [1], the entanglement of a finite region of size L 



with the remainder of the system grows logarithmically 
in L at a critical point, while away from criticality the 
entanglement is localized near the boundaries of the sub- 
system and goes to a constant for large L. For critical lat- 
tice models like quantum spin chains, the entanglement 
of a segment of L sites with the remaining sites grows 
as log2 L, with a coefficient determined by the central 
charge of the conformal field theory (CFT) [5] : 
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■ log2 L. 



(3) 



Here c and c arc the holomorphic and antiholomorphic 
central charges of the CFT (for the cases we discuss 
c — c), which control several physical properties such 
as low-temperature specific heat. Slightly off criticality, 
the spin chain has a finite entanglement length ^, and the 
entanglement saturates as i — > oo to 5 ~ ^^ log2 S,. 

An example of a quantum spin chain that is criti- 
cal without disorder is the antiferromagnetic Heisenberg 
model; the ground state of the spin-half chain 



i/ = J^S,.S,+i = J^(5f5| 



s+s:: 



-^^) (4) 



is quantum critical for the antiferromagnetic case J > 0. 
Staggered spin-spin correlations, (—1)1*^^1(8^ • Sj), fall 
off as l/|z — j'l up to logarithmic corrections. 

The nature of quantum spin chains with quenched ran- 
domness at zero temperature is quite different from the 
above pure case. It is believed that any initial random- 
ness in the distribution of couplings drives the system at 
long distances to a random quantum critical point: in 
RG language, disorder is a relevant perturbation to the 
pure critical points. This flow to strong disorder occurs 
for the Heisenberg chain, the XX chain, which has cou- 
pling only in two spin directions, and the quantum Ising 
chain, which has couplings in one spin direction plus a 
normal magnetic field (both made random). 

The low energy properties of the random Heisenberg 
and XX models are described by the random singlet 
phase [15]. This is shown using the real-space RG ap- 
proach [3, 9, 10]. We review the real-space RG approach 
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FIG. 1: a. RG decimation step: if Ji is the strongest bond in 
the chain, sites i and i + 1 form a singlet (solid line); which 
diagonalizes JiSi ■ Si+i. Quantum fluctuations produce an 
effective interaction between sites i — 1 and i + 2 (dashed 
line), b. The random singlet ground state. Note that singlets 
may connect arbitrarily distant sites. 



and the random singlet ground state, starting with the 
random Heisenberg Hamiltonian 



H 



J2 J^^^ 



'i+l- 



(5) 



The same results apply to the XX chain. In Eq. (5), J^'s 
are drawn from any nonsingular distribution [3] . 

The real-space RG analysis consists of iteratively find- 
ing the strongest bond, e.g. Ji, and diagonalizing it inde- 
pendently of the rest of the chain. This leads to a singlet 
between spins i and i -t- 1 in zeroth order (Fig la): 

|V(°)) = IV'.<«>-^ (IT4m) - U.Wi) ) IV^.>m)- (6) 

Next, we treat the rest of the Hamiltonian as a perturba- 
tion. If we begin with strong disorder (the distribution of 
In Ji is wide), we can assume that Ji ^ Ji-i, Ji+i, and 
use degenerate second order perturbation theory. This 
leads to a Heisenberg interaction between the neighbor- 
ing spins at sites i — I and i + 2 with strength 



J, 



i-l,i+2 



Jj-lJj- 

2Ji 



< Jt- 



Ji, Ji+1- 



(7) 



Thus we eliminate two sites, and reduce the Hamilto- 
nian's energy scale. Iterating these steps gives the ground 
state. Although this method is patently not correct when 
applied to a chain with little disorder, it is still applicable 
and is asymptotically correct at large distances. [3] 

The RG leads to an integro-differential flow equation 
for the bond coupling distribution. This equation is 
best stated in terms of the logarithmic coupling strength 
(3 = \nj and RG flow parameter F = In ^ , where Uq 
is the Hamiltonian's initial energy scale, and O is its re- 
duced energy scale. These variables capture the scaling 
properties of the problem; e.g., neglecting a In 2, Eq. (7) 
is simply /3i-i,i+2 = A-i + Pi+i- Note that strongest 
bonds have (3 = 0. In terms of /3 and F we have [3] 
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(8) 



FIG. 2: The entanglement entropy of a segment is the number 
of singlets that connect the segment with the rest of the chain 
(shaded area). In this example there are two such singlets. 



The following solution is an attractor to essentially all 
initial bond distributions, and it describes the low-energy 
behavior of the spin chain: [3] 
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1 



-/3/r 



(9) 



This is the random-singlet fixed point distribution. 

The real space RG shows that the spin chain is in the 
random-singlet phase. In this phase singlets form in a 
random fashion over all length scales, and can connect 
spins arbitrarily far apart (Fig. lb). Long-distance sin- 
glets form at low energy scales. On average, 



A 



(10) 



where A is the length scale of singlets forming at energy 
scale il = rioE^^ . The long range of the low-energy sin- 
glets leads to slowly decaying average correlations, which 
for the random Heisenberg model decay algebraically, 
and not exponentially as one would expect from the lo- 
calized nature of the random-singlet state. 

Let us focus on the entanglement entropy of the ran- 
dom Heisenberg model. The entanglement of a spin-1/2 
particle in a singlet with another such particle is 1 , which 
is the entropy of the two states of a spin with its partner 
traced out. The entanglement of a segment of the ran- 
dom Heisenberg chain is just the number of singlets that 
connect sites inside to sites outside the segment (Fig. 2). 

To obtain the entanglement, we calculate the number, 
N, of singlets that form over a single bond B (in the 
example above, we form singlets between sites i and i + 
1, and later in the RG between sites i — 1 and i + 2, 
etc.). If (as a first approximation) we neglect the history 
dependence of the distribution of bond B, we can find 
N by using the distribution of bond strengths, Eq. (9). 
When we change the energy scale ft ~* ft — dil, F -^ 
F -I- dF, all bonds with < ^9 < dF get decimated. The 
average number of decimations over bond B grows by 



7T^ 

dN = dTP{(3 = 0) = -fT 



(11) 



which leads to iV = InF. 

This picture breaks down once singlets exceed the size 
L of the segment , when (by Eq. 10) F = \/X. So the 
number of singlets emanating from the segment of size 
L, i.e. its entanglement with the rest of the chain, is 



S'l ~ ^L « 2 • In VT + fc = In L + A:, 



(12) 



where A; is a nonuniversal constant, which also depends 
on the initial realization of the disorder. 

Neglecting the history of B allowed us to see simply 
why the entropy depends on InL, but the coefficient we 
obtained is not correct. To include the history of i3, we 
note from Eqs. (7, 9) that the bond strength distribution 
of B right after being decimated at Fq is 



Q(/3) 



dMf32Sif3i+p2-0)Pro {P2)Pro [Pi 



"rf ■ 

(13) 

Now we ask at what F is B decimated again. To answer 
this, we need a flow equation for Qr{P) similar to Eq. (8). 

In the following, we use the convention that / d/3Qr (/3) — 



pr is the probability that bond B was not yet decimated 
at scale F. With this convention in mind, Qr{P) obeys 








u^-<|>H^ 


^=B 


' 


■^ 



FIG. 3: (a) Typical ground state in the random quantum 
Ising model. It formed as follows; sites 2 and 3 form a cluster, 
site 1 is frozen in the x direction, site joins the cluster of 
2 and 3, and finally the large cluster is frozen along x. (b) 
The entanglement of a segment L is given by the number of 
decimated clusters that connect the segment with the rest of 
the chain (shaded area). In this example there are two such 
clusters: sites 2 and 3, and sites 1 and 4. 
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2Qr(/3)Pr(0) 



^2Fr(0) / rf/3id/32<5(/3i + /32 - /?)Pr(/32)Qr(/3i 



(14) 



The first term is due to the change in /3 when fJ changes, 
the second and third terms account for i?'s flow due to 
one of its two neighbors forming a singlet. Note that 
-^ = — Qr(0). Eq. (14) can be solved using the ansatz: 



We discuss the interpretation and physical conse- 
quences of this effective central charge below, but first 
obtain its value for the quantum Ising case. The pure 
quantum Ising model has a well known ferromagnetic to 
paramagnetic phase transition. Its random analog is 



H 



^ Jia^cFl^i 



^/l,<Tf, 



(19) 



QAP) 



av + hr^]Pr{p) 



(15) 



by substituting Eq. (15) in Eq. (14) we obtain 
i^ = F^ = 6r - 2ar, %^ = T^ = -6r 



di 



dV 



ar, 



(16) 
with I = InF/Fo. Also apo = 0, bpo = 1, from Eq. (13). 

Next we calculate the rate of singlet formation over 
B. First, note that the survival probability pr obeys 
pr = Or + br and depends on F only through I = InF/Fp. 
Therefore the duration / between two consecutive sin- 
glets forming on bond B is independent of Fq. This also 
proves that the number of singlets over B is proportional 
to InF. To find the proportionality constant we calcu- 
late the average duration {I) between decimations; the 
number of bonds is then N = ^j^. We have 



oo 

{I) = / dpr I ~ / dlarl- 



(17) 



From Eq. (16) one finds ar 
Inserting this in Eq. (17) we find (1) 



1 



3-v/5 , 



3+v/5 , 



SL = ^-2\nT + k=^-^log^L 



3. Therefore: 

k. (18) 



Hence the 'effective central charge' of the random Heisen- 
berg chain is c = 1 • ln2, which is the central charge of 
the pure Heisenberg chain times an irrational number. 



where hi and Ji are random, and a are spin-1/2 Pauli 
matrices. This model also has a ferromagnetic to param- 
agnetic phase transition, described by a random critical 
point similar to the random singlet phase [4]. Hence we 
expect entanglement in the quantum Ising case also to 
scale logarithmically with the size L of the test segment. 
As with the random Heisenberg case, we use real-space 
RG to study the ground state of the random quantum 
Ising model (Fig. 3a). Again we diagonalize the term 
with the largest energy scale in the Hamiltonian; if it is 
~Ji<7iaf_^_i, we set sites i and i -I- 1 to point in the same 
direction, \ip)i,i+i = a |TiTi+i) +/? |iiii+i), thus creating 
a ferromagnetic cluster. Quantum fluctuations yield an 
effective transverse field on the cluster, 

hij+i = -^ — (< hi, hi+i, Ji) . (20) 

If the term with largest energy happens to be —hiaf, 
we set the i'th spin point in the x direction, 1-0)^ = 
|— ^j) — -y= (Ifj) -I- lli)), by which we decimate the I'th 
spin. Quantum fluctuations produce an effective Ising 
coupling between sites i — 1 and i + 1 with strength: 

J^-l,^+l = ^^^ (« J.-l, Jz, h,) . (21) 

The RG flow equations for the distributions of hi and 
Ji as a function of fi = max{ J,;, hi} support an attractor 
in which the logarithmic coupling distributions, i?(C) and 
P(/3), withC = ln^, /? = 
singlet expression in Eq. 



In 



n 



■^, are given by the random 
(9). As for the random singlet 



scaling, at criticality the absolute length A of the domains 
decimated by the transverse field scales as A ~ F^. 

At low energy, larger and larger ferromagnetic clus- 
ters are formed and then decimated. For a segment of 
length L, ferromagnetic clusters which are completely 
within or completely outside the segment and decimated 
by the transverse field do not affect the entanglement. 
The only contributions come from ferromagnetic domains 
that cross the boundary of the segment (Fig. 3b), and 
each such cluster contributes 1 to the entanglement. 

Next we calculate how many ferromagnetic clusters are 
formed and decimated over an edge of the segment. At 
a given energy scale, the edge of the segment can either 
separate two unpaired sites, or be contained in a cluster 
which is partially in the segment (when such a cluster 
is decimated the edge returns to the first case). At the 
critical point these possibilities must occur with the same 
probability by a self-duality of the quantum Ising model 
([4, 11]). Hence the probability that the edge is in an 
active cluster at scale F is p = 1/2. By the analysis as 
for the Heisenberg chain, we obtain A^(r) = p- 5 InF, and 



Sl 



■InL 



In 2 



■ log2 L 



(22) 



with k a non-universal constant. The effective central 
charge of the random quantum Ising model is c = 1/2-ln 2 
- In 2 times the central charge in the pure system. 

Eq. (22) shows that the critical quantum Ising chain 
has half the entanglement of the random singlet phase. 
But both of these are in the same infinite-randomness 
fixed point scaling category, so a difference in the entan- 
glement entropy may seem surprising. However, these 
two systems also differ in their temporal correlations: 



< a^{0)a^„{T) > 



1 



In^ OqT 



(23) 



with z = 2 for the XXZ models and z = 1 for the quan- 
tum Ising model. These two similar strong-randomness 
fixed points belong to different universality classes, and 
the effective central charge is sensitive to this difference. 

For pure chains, the prefactors of correlation functions 
are non-universal but the central charge is universal. In 
the random case, non-universal correlation prefactors are 
generated by inaccuracies of the order of the lattice spac- 
ing in the location of the low energy effective spins; these 
occur when the RG scale fl is still large [12]. Such errors 
do not affect the universal coefficient of the logarithmic 
divergence of the entropy, but modify the additive non- 
universal 'surface term' k in Eqs. (18, 12). Note that 
entanglement is self- averaging as chain length N ^ 00. 

The central charges we find for the random Heisen- 
berg, XX, and quantum Ising chains are those of the 
pure models times In 2. Although irrational, these cen- 
tral charges are universal quantities which describe the 
universality class of the random chains. An example of 
the importance of c in the pure case is the well known 



"c-theorem" [13] that if an RG flow connects one criti- 
cal point A to another critical point B, then ca > cb- 
The values of c obtained here for random systems sug- 
gest, given the relevance of disorder in these systems, 
that there is a generalized c-theorem based on entangle- 
ment even for nonconformal quantum critical points in 
ID. This may imply constraints on the values of 2 for spin 
chains obtained by, e.g., disordering higher-spin CFTs. 

Since only certain rational values of c are allowed for 
well-behaved CFTs with c < 1, the irrational c for ran- 
dom critical points is a fundamental difference between 
pure and random cases. This is reminiscent of the irra- 
tional residual entropy that appears in quantum impurity 
problems and satisfies a "(/-theorem." [14] 

The ratio between the random and pure values of c 
is unexpectedly the same for all the different chains we 
studied. Perhaps c of any random fixed point derived 
from a pure conformally invariant point is the product 
of the central charge of the pure theory and a universal 
number determined by the flow from the pure to the ran- 
dom fixed point. Numerics on the random XX model via 
its free-fermion representation could determine whether 
c appears as a universal amplitude as in the clean case 
for quantities beyond entanglement. It is clear already 
that the universal logarithmic scaling of entanglement 
provides a powerful way to characterize both pure and 
random quantum critical points in one dimension. 
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